Probing singularities in quantum cosmology with curvature scalars 
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We provide further evidence tliat tlie canonical quantization of cosmological models eliminates 
the classical Big Bang singularity, using the DeBroglie-Bohm interpretation of quantum mechanics. 
The usual criterion for absence of the Big Bang singularity in Friedmann-Robertson- Walker (FRW) 
quantum cosmological models is the non-vanishing of the expectation value of the scale factor. We 
compute the 'local expectation value' of the Ricci and Kretschmann scalars, for some quantum FRW 
models. We show that they are finite for all time. Since these scalars are elements of general scalar 
polynomials in the metric and the Riemann tensor, this result indicates that, for the quantum 
models treated here, the 'local expectation value' of these general scalar polynomials should be 
finite everywhere. Therefore, according to the classification introduced in Refs. we have 

further evidence that the quantization of the models treated here eliminates the classical Big Bang 
singularity. 
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The presence of singularities in cosmological models in general relativity is an old issue 11]. In a series of articles. 
Hawking, Penrose and Geroch showed that, if certain very general conditions are satisfied, singularities are always 
present in cosmological models based on general relativity [1^ . One of such singularities is the Big Bang singularity, 
which is believed to represent the very beginning of the Universe. Here, one has a fundamental problem because, if the 
beginning of the Universe is a singular event of general relativity, that theory cannot describe it. In order to overcome 
that fundamental problem, many authors proposed the quantization of gravity. Quantum cosmology was the first of 
such attempts and, since the first model, has showed good signs toward the solution of the above mentioned problem 
[l|. Since then, many important works have been done by computing the wave function of the universe (4') in different 
minisuperspace models. Some authors find 5* by solving the Wheeler-DeWitt equation 5, 6, 13 18, 24] and others by 
using the path integral approach (T9l - |2T| . As a common result, in most of them the problem of the initial singularity 
was claimed to be solved. The main argument used to support those claims depends on the quantum mechanical 
interpretation used in each particular work. 

The two interpretations most frequently used in quantum cosmology are the Many Worlds one [2^ and the 
DeBroglie-Bohm one 0- As in the usual Copenhaguen interpretation of quantum mechanics, in the Many Worlds 
interpretation one cannot talk about trajectories of the canonical variables, but only about mean values of those 
variables. On the other hand, in the DeBroglie-Bohm interpretation the trajectories of the canonical variables are 
meaningful and can, in principle, be computed by solving a system of differential equations involving derivatives of the 
wave function phase. In the minisuperspace models [3] treated using the Many Worlds interpretation, the common 
argument used to justify the absence of a Big Bang singularity is the fact that the mean value of the scale factor 
(a), as a function of a chosen time, never vanishes 0, [13, [3 In the models investigated according to the 
DeBroglie-Bohm interpretation, the argument was that the scale factor Bohmian trajectories a(t) as a function of a 
chosen time never go through a = [^-Q. That result is supported by the fact that the quantum potential present 
in the dynamical equation of a, for those models, is repulsive for a near to zero 043- In the present work, we shall 
restrict our attention to the DeBroglie-Bohm interpretation of quantum mechanics. 

Quantum cosmology was the first attempt in order to remove the Big Bang singularity by quantizing the gravita- 
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tional theory [ij. The DeBroglie-Bohm interpretation of quantum mechanics [3|, is frequently used in quantum 
cosmology. In the minisuperspace models treated using the DeBroglie-Bohm interpretation the common argument 
used to justify the absence of a Big Bang singularity is the fact that the scale factor Bohmian trajectories a{t), as a 
function of a chosen time, never go through a — 043 ■ 

In order to derive the physical content of any operator A(x,px), in the DeBroglie-Bohm interpretation of quantum 
mechanics, one must compute the so-called 'local expectation value' of that operator defined as [1] 



'i'*{x,t)-^{x,t) ] ' ^ ' 

where {A'i'){x,t) = J A{x,x')'^{x' ,t)(fx' . If we apply this definition in quantum cosmology, it is easy to see that 
the 'local expectation value' of the scale factor operator is the real, time dependent, scale factor function. On the 
other hand, the 'local expectation value' of more complicated operators constructed out of the scale factor and its 
canonically conjugate momentum will be much more difficult to compute and in general will require a specific factor 
ordering prescription. 

The Big Bang singularities that occur in the classical Friedmann-Robertson- Walker (FRW) models are said to be 
'scalar polynomial singularities' [9]. A 'scalar polynomial singularity' is the end point of at least one curve on which 
a scalar polynomial in the metric and the Riemann tensor becomes infinite [loj . In the present work, we wish to give 
further evidence, besides the usual one, that the canonical quantization of FRW models removes the initial Big Bang 
singularities of those models. We shall compute, for all FRW models considered, the 'local expectation value' of the 
Ricci (= ga/sg-ysR^"'^^) and the Kretschmann scalars (= Rap-ysR"'^'^^)- As we shall see, they are finite for all time. 
Therefore, since these scalars are components of general scalar polynomials in the metric and the Riemann tensor, 
that result indicates that the 'local expectation value' of a general scalar polynomial in the metric and the Riemann 
tensor should be free from singularities in those models. 

In the present work, we shall consider FRW cosmological models coupled to a radiative perfect fluid, treated by 
means of the variational formalism developed by Schutz [23l] . Our main motivation to choose the matter content of 
the model as radiation is because we would like to describe the very early Universe, the so called 'radiation dominated 
era'. At that time, the quantum effects were more important. The models are described by the Hamiltonian @ 



A{x,t) = Re 



n = !^+6ka'-PT, (2) 

where pa and pt are, respectively, the momenta canonically conjugate to the scale factor (a) and the radiation variable 
(T). The parameter k is related to the spatial curvature of the model and may assume the values +1 (positive 
curvature), —1 (negative curvature) and zero (no curvature). We employ the natural system of units, defined by 
h = c= IGttG = 1. 

The quantization of those models follows the Dirac formalism for quantizing constrained systems. The application 
of this formalism for the present models result in the following Wheeler-DeWitt equation for the wave function \l/(a, T) 

0'. . 

First we introduce a wave function which is a function of the canonical variables a and T, 



^ = *(a,r) . (3) 

Then, we impose the appropriate commutators between the operators a and T and their conjugate momenta pa 
and Pt- Working in the Schrodinger picture, the operators a and T are simply multiplication operators, while their 
conjugate momenta are represented by the differential operators 

d d 

Pa ^ , PT ^ . (4) 

Finally, we demand that the operator corresponding to T-L annihilate the wave function 5*, which leads to the 
Wheeler-DeWitt equation. From Eq. ([2]), this formalism leads to 

— - 144fca2$ + 2At— = 0, (5) 
oa'' OT 

where T = — t. Several solutions to equation ([5]) are known. 
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Consider first the case k = 0. Then the solution to Eq. ([5]) is given by [23 |. 



4'(a, r) 



4cr 



1/4 



6i 

(2crr-z)(2r + i) ^ I r 



1 + 



2ctt — « 

In the cases k — ±1, the solution of Eq. ([5]) may be written in the following form Q: 



(6) 



X exp ■ 



1/4 



A: 



cos^ (Vkr) [2a tan(\/fcT) — [2 tan( VXt) + i^/k] 



6iVk 
tan(-\/A:T 



iVk 



cos^ {\/kT)[2a ta,n{^/kT) — i^/k] 



In order to use the DeBroglie-Bohm interpretation we must rewrite ^' in the polar form 



*(a,r) 



4cr 



1/4 



e exp (iS) 



Consider first the case k = 0. Then the solution to Eq. ([5|) is given in polar form by [2J], 



(7) 



(8) 



e = goir) X exp 



-12cra^ 



1 + 4ct2t2 



and 



(9) 



S = Mr) + 



24crV 



For the cases fc = ±1 one has 



(10) 



= 9k{T) X exp 



-Uaka'^ 



k cos^ (Vkr) + 4cr2 sin^ (Vkr) 



and 



(11) 



S = /fc(r) + 



2tan(yfcT) 



1 - 



k cos^ (Vkr) + 4a-2 sin^ (%/A:r) 



(12) 



Following the DeBroglie-Bohm interpretation and the fact that Pa = 12d, where the dot means differentiation with 
respect to r, we may compute the scale factor Bohmian trajectory from Q 



Pa 



dS_ 
da 



(13) 



For the case fc = 0, from Eqs. ® and (HU]) we have [1, H] 



a(T) =ao(l + 4aV^) 



2 2n1/2 



where ao is an integration constant and 



Q 



l + 4cr2r2 (l + 4cr2r2)2 

For the cases k = ±1, Eqs. ^ and ^ yield 



(14) 



(15) 



a{T) = ak{k cos^ (Vkr) + Acr^ sin^ (Vkr))^/^, 



(16) 
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where is an arbitrary integration constant, associated to a given k, and 

ak 24cr2fc2a2 



k cos^ (\/fcT) + 4ct2 sin^ (Vfcr) (fc cos^ {^/kr) + 4ct2 sin^ (VfcT))^ 



(17) 



In all these cases it is clear that the scale factor Bohmian trajectories never reach a = 0. Therefore, one may 
conclude that these models are free from the Big Bang singularity. 

Here, we would like to give further evidence, besides the non-vanishing of the scale factor, that those models are 
free from the Big Bang singularity at the quantum level. In order to do that, we shall compute the 'local expectation 
value' of the Ricci and the Kretschmann scalars. As it will be seen, they remain finite for all time. These scalars are 
components of general scalar polynomials in the metric and the Riemann tensor. Therefore, we have an additional 
indication that the 'local expectation value' of these general scalar polynomials should be finite everywhere. 

For the present models, the expressions of the Ricci (i?) and the Kretschmann {K) scalars are given, respectively, 
by 

R = g''^g'"'Rai,cd = ^ + 5 and (18) 



a2 



7^ jyabcdry Pa PaPa , Pt , ^Pa , 12fc^ 

where we used Pa = 12d in order to write R and K in terms of Pa, the momentum canonically conjugate to a. 

It is important to notice, before we proceed, that R and K given by Eqs. (|18l) and ([T51) are to be promoted to 
quantum operators. Since, at the quantum level, Pa and a do not commute, we shall have to introduce a specific 
factor ordering in order to correctly describe the terms involving products of powers of a and Pa ■ Here, we shall use 
a symmetrization procedure known as the Weyl ordering [25J. In order to obtain the Weyl-ordered expression of a 
product (a"p™), one first randomly orders the a's and p'^^s, with each different ordering counted once, then divides 
the result by the number of terms present in the final expression (25j . 

In Eqs. ()18p and (|19|) we notice the presence of the time rate of change of the momentum pa- Quantum mechanically, 
it is an operator and in the DeBroglie-Bohm interpretation it has the following value [sj, 

Pa^-y{V + Q) (20) 

Where V is the classical potential present in the Hamiltonian ^ and Q is the quantum potential. From Eq. ([2]) V 
is given by 

V = 6ka^ . (21) 
Using expressions (fTSl) and p7)) for Q, and Eq. ((2T|| for V, we may compute the time rate of change of the momentum 

Pa- 

For fc = we obtain 

(22) 

while for fc = ±1 we find 



Pa^~12ka+ (23) 

Introducing the above values of Pa in the expressions for the Ricci and Kretschmann scalars we obtain new forms 
which depend only on the operators a and Pa- 

For fc = 0, introducing pa given by Eq. ([^ in Eqs. and ([TO)), we get 

« = ,24, 

al (1 + 40-2x2)3' 
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1 pi 2aia^ a x ^ , . 



For k = ±1, inserting pa from Eq. ([23]) into Eqs. (fT8| and p9|) . we obtain 

i? - ^^"^^ f26) 
(fc cos2 (v^r) + 4a2 sin^ {VkT))^ ' 

K=^4- ^4 + ^4 + 192a V,4i - 96a2fc3at4 + 24^ . (27) 
864 3 3 "a^^ ^ ^ 

The expressions for R both in the cases k = and k = ±1 depend only on the operator a. Therefore, the R 'focal 
expectation value', for both cases, is given by the same expression of the operator R with the operator a replaced 
by the real, time dependent scalar factor function. It is clear from the above expressions for R that it is regular at 
T = 0. Thus, R is regular for all values of k at the beginning moment of the universes described by the corresponding 
models. In order to compare the behavior of the R 'local expectation value', as a function of r, with the classical 
expression of R, we produced Figs. ([T]), ([3|) and (O, one for each value of k. The classical scalar factor was derived 
with initial conditions compatible with those of the scale factor Bohmian trajectories. One may easily see from those 
figures that, for the cases where the spatial sections are open, both quantities coincide for large t. 

In order to derive the physical content of the operator K, in the DeBroglie-Bohm interpretation of quantum 
mechanics, we must compute its 'local expectation value' Eq. ((!]). Due to the presence of the terms 



4 2 2 

Pa Pa Pa 



(28) 



in the expressions of K, Eqs. (|25p and (P?]). we choose to use the Weyl ordering '29]. 

In order to reduce the Weyl-ordered expressions of each one of the above products of a's and p^s, we use the 
commutation relation between the operators a and Pa- With its aid, we write all terms, in each of the Weyl-ordered 
expressions, such that pa or a power of Pa must appear to the right of a or a power of a. This procedure simplifies 
very much the Weyl-ordered expressions because most of the terms combine with each other. With the aid of the 
following commutators, 

[o",p^] = 2ina"~^pa + n{n-l)a"-'^ 

[a",p^] = Aina"~^pl + 6n{n-l)a''-^pl~4in{n-l){n-2)a"-^Pa 

- n{n-l){n~2){n~3)a''-^, (29) 

where n is a positive or negative integer, we obtain the following Weyl-ordered expressions for each one of the products 
of a's and p'^s in Eq. ((28t : 



,10 



W 



,6 



W 



76^'a + 6^-yPa- 23 — . (30) 



Now, we substitute Pa given by —id/da into these Weyl-ordered expressions and compute their 'local expectation 
value' Eq. (P), using the wave function (jSj. It is important to remember that the wave function ^ must be written 
in the polar form ([SJ. Then, we obtain the expressions for the 'local expectation values' of each operator in Eq. ([50)1 
as functions of 6(a, r) and derivatives of Q{a, r) and 5(a, r) with respect to a. 

pi\ _ {dS{a,T)/da)^ ■i{d'^S{a,T)/da^Y ll&{dS{a,T) / daf 



6 



+ 



nW 



12(99(a, T)/da){d^S{a, T)/da^)(dS{a, T)/da) 

a8e(a,T) 
48(5^ 5'(a, T)/da^){dS{a, T)/da) 

&{d^Q{a,T)/da^){dS{a,T)/daf _ l&{d^Q{a,T) / da^) 

a8e(a,T) a9e(a,T) 
116(9^e(a,T)/9a^) 9-^9(a, r)/9a4) _ 440(99(a, T)/9a) 

aioe(a,T) ^ a8e(a,T) a^^Q{a,T) 
48(99(o, r) /9a)(c)5'(a, T)/Oa)^ _ Ajd"^ S{a,T) / da'^){dS{a,T)da) 
a9e(a,r) ^8 

10831 
15^' 

175 (a5'(a,r)/aa)2 



10(ae(a,r)/aa) (929(0, r)/aa2) 



23 
2^ 



6(9e(a, r)/9a) 
0^9(0, r) 



a"9(a,T) 

(9^9(0, T)/9a2) 
a69(a,T) 



ai09(a,r) 
(aS'(a,r)/aa)2 



(31) 
(32) 

(33) 



Finally, we compute the K 'local expectation value' for both cases of fc = Eq. ((25t and k = ±1 Eq. p7|) . Then, 
in order to compute the K 'local expectation value', After that, we introduce the values of 9(a. r) and S'(a, r) for 
each k, Eqs. ([MT^ in the 'local expectation values' of {p^/a^)w: {Pa/a^^)w and {p'^/a^)w 

Finally, we combine them following Eq. ([25]) for k = and Eq. (P7)) for A; = ±1. 

For k = we find 



where 



For fc = 1 we get 



,12 



(l+4CT2r2)^ 



6144aV'* - 10752aV2 + 960cr'*, 
-1920crV + 304cr3, 



4^2 



-56crV 

9cr, 

10831 



601 
~9~ 



12960 



(34) 



(35) 



where 



a];2 (cos2 r + 40-2 sin2 t)^ 



24(256cr^ + 256cr^ - 160cr'* + 16cr2 + 1) cos^ t 

96cr2(-128cr^ - 16cr'* + 40cr2 - 7) COs2 T 
192cr'*(32cr4-24cr2 + 5), 

80cr(16cr'' - 8cr2 ^^^4 ^ _^ 40cr(-16cr'* + 16cr2 

3)cos2r- 16cr^(40cr2 + 19), 

(17/15)(-4(t4 + 2^2 - (1/4)) cos^ r + (l/5)((988/3)fT4 
(34/3)cr2 - (71/4)) cos2 T + (2/3)f72(_92cr2 + (757/15)), 
(27/10)a, 
(10831/12960). 



(36) 



(37) 
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Finally, for = — 1 we obtain 



a^^-^ (cosh^ r + 4(t2 sinh^ t)^ 



(38) 



where 




— 1 



24(256cr** - 256cr^ - IGOcr'' - 16cr^ + 1) cosh"^ r 
96(T2(-128cr'^ + IGcr^ + 400-2 + 7) cosh^ r 
192cr^(32cr* + 24cr2 + 5), 

80cr(-16cr'' - SfJ^ - 1) cosh^ T + 40cr(16cr'' + IGcr^ 
3) cosh^ T + 16cr^(40cr2 + 19), 

(17/15)(4cr^ + 2^2 + (1/4)) cosh* r + (1/5)(-(988/3)(t* 
(34/3)cr2 + (71/4)) cosh^ r + (2/3)cr2(92cr2 + (757/15)), 
(27/10)0, 
(10831/12960). 



+ 



+ 




1 



+ 



1 




1 



1 



(39) 



It is clear from Eqs. (IM)) . ([55]) and ([551) that the A' 'local expectation value', for each model, is regular for all r. 
Including the limit t — > 0, that is, at the beginning moments of the corresponding classical universes. In order to 
compare the behavior of the K 'local expectation value', as a function of r, with the classical expression of K , we 
produced Figs. ([2]), @ and ([S]), one for each value of k. The classical scalar factor was derived with initial conditions 
compatible with those of the scale factor Bohmian trajectories. One may easily see from those figures that, for the 
cases where the spatial sections are open, both quantities coincide for large t. It is important to notice, that this 
result is independent of the factor ordering used here. In fact, observing Eqs. (pi]) . ([36]) and ([38]) . we conclude that 
they are regular mainly because the scale factor Bohmian trajectories a{t) never go through a = 0. Therefore, from 
the operatorial expression of K ([TT)) . it is not difficult to see that whatever factor ordering we decide to use the 
denominator of the K 'local expectation value' will be a polynomial in the scale factor. Then, if we take in account 
that it does not vanish for any t, the K 'local expectation value' will always be regular. Now, since, R and K are 
elements of general scalar polynomials in the metric and the Riemann tensor, the above results indicate that, for 
the quantum models treated here, the 'local expectation value' of these general scalar polynomials should be free of 
singularities. Therefore, according to the classification introduced in Refs. [1, [13]) we have further evidence that the 
quantization of the models treated here eliminates the classical Big Bang singularity. 

We believe that the above result may be extended to FRW models with matter contents described by other types 
of perfect ffuids. This is the case because, as we have mentioned above, the main reason for the regularity of the 
'local expectation values' of R and K is that the scale factor Bohmian trajectories a(t), as a function of a chosen time, 
never go through a = 0. Therefore, if we consider FRW models with matter contents described by other types of 
perfect ffuids and we obtain Bohmian trajectories a{t) that never go through a = 0, very likely the 'local expectation 
values' of R and K of those models will always be regular. In Reference (5j], the authors calculated the scale factor 
Bohmian trajectory, as a function of the proper time, for a flat (fc = 0) FRW model with matter content described 
by a generic perfect fluid. A generic perfect fluid is described by the equation of state p = wp, where w is a constant. 
In the radiation case w = 1/3. There, they found that the scale factor never go through a = whatever value w 
assumes. Therefore, for that case we are very confldent to say that the 'local expectation values' of R and K should 
always be regular. Unfortunately, we do not have such general result for models with k = ±1. On the other hand, if 
one believes that the quantization of those models will solve the singularity problem, they will also have scale factor 
Bohmian trajectories that never go through a = 0. In this way, the 'local expectation values' of R and K of those 
models should, also, always be regular. 
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FIG. 1: The R 'local expectation value' for fc = 0, cr = 1 
and ao = 1. It is always regular and goes to zero for large 
r. For this case R evaluated over the classical scale factor is 
identically zero. 



FIG. 2: The K 'local expectation value' for k — Q, a = 1 
and oo = 1. It is always regular and goes to zero for large 
T. For this case K evaluated over the classical scale factor is 
identically zero. 
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FIG. 3: The lower curve represents the R 'local expectation 
value' for k — 1, a — 1 and ai = 1. It is periodic and al- 
ways regular. The upper curve represents R evaluated over 
the classical scale factor. 



FIG. 4: The lower curve represents the K 'local expectation 
value' for = 1, a = 1 and ai = 1. It is periodic and al- 
ways regular. The upper curve represents K evaluated over 
the classical scale factor. 




FIG. 5: The upper curve represents the R 'local expectation 
value' for k = —1, a = 1 and a-i = 1. It is always regular. 
The lower curve represents R evaluated over the classical scale 
factor. Both curves go to zero for large r. 



FIG. 6: The lower curve represents the K 'local expectation 
value' for k = —1, a = 1 and a_ = 1. It is always regular. 
The upper curve represents K evaluated over the classical scale 
factor. Both curves go to zero for large t. 



